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The kinetic theory of evaporative cooling developed by Luiten et al. [Phys. Rev. A 53, 381 (1996)] 
' is extended to include the dependence of the elastic scattering cross section on collision energy. We 

introduce a simple approximation by which the transition range between the low-temperature limit 
and the unitarity limit is described as well. Applying the modified theory to our measurements on 
evaporative cooling of metastable helium we find a scattering length \a\ = 10(5) nm. 
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0\ ', I. INTRODUCTION 



Evaporative cooling is at present the most powerful method to increase the phase-space density of a dilute trapped 
gas in order to reach Bose-Einstein condensation (BEC) or Fermi degeneracy 0. The atoms are first laser-cooled 
to temperatures of 1 mK or lower and trapped either magnetically in a Ioffe quadrupole (IQ) trap, or optically in a 
dipole trap. Here the energy of the atoms is redistributed by elastic collisions. A few atoms acquiring energies larger 

S than the trap depth e t (with a maximum energy of 2e t ) are expelled from the trap. This reduces significantly the 
mean energy of the trapped ensemble while only a relatively small number of atoms is lost. In forced evaporative 
"■^J , cooling the trap depth is gradually lowered, during which cooling increases the density and the elastic collision rate, 
£h ' accelerating evaporation. 

In the Walraven group kinetic theory has been developed to describe this evaporative cooling process for cold atomic 
, i eases 00. The theory assumes collisions occurring in the s-wave limit, an energy-independent clastic scattering cross 
section a and sufficient ergodicity in the trap. A truncated Boltzmann distribution is introduced and shown to be 
t-H , consistent with the classical Boltzmann equation applied to a trap of finite depth. Then equations of statistical 
mechanics are obtained describing the kinetics, particle loss and energy loss. Closed expressions are derived that are 
directly applicable to an IQ trap. However, the assumption of an energy-independent cross section a is not justified 
when evaporation starts in the transition range between the low-temperature limit and the unitarity limit. This 
\ is the case, for instance, in metastable helium, where evaporative cooling typically starts at ~1 mK 0, Hj,El- I n 
■ order to understand evaporative cooling under these conditions as well, this paper extends the theory to incorporate 
an energy-dependent a. The theory is subsequently applied to interpret our evaporative cooling experiments with 
metastable helium, allowing the s-wave scattering length a to be extracted. As BEC has been realized by two groups 
in France Q, 0, a is known to be large and positive. Published experimental values are a = +20(10) nm Q and 
a = +16(8) nm both obtained from Bose-condensed clouds. Recent experiments in Paris .7] with clouds just 
above the BEC transition suggest that a is smaller than previously published values. Theoretical values determined 
from precisely calculated molecular potentials are +8.3 nm Q and +12.2 nm |9LIToj. 

This paper is organized as follows. In Section[H]the theory of evaporative cooling based on the papers of Luiten et 
al. and Pinkse et al. is summarized and extended. First, Section III Al gives the density and energy of an atomic 
cloud in a magnetic trap of finite depth and Section III Bl provides the equations necessary to calculate all relevant 
loss rates during the cooling process, assuming a is energy- independent. Then, in Section III CI the atom and energy 
loss rates due to evaporation are recalculated, now with an energy-dependent cross section a. Additionally, a simple 
approximate expression is presented for an effective cross section that can be implemented directly in the equations 
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for the evaporation loss rates given before in Section III Bl The theory is used in Section II I II to interpret evaporative 
cooling experiments in metastable helium performed in Amsterdam and a value for the scattering length is deduced. 
Section ll V I summarizes and compares the extracted scattering length with results obtained in the other experiments 
and with theoretical values from molecular potential calculations. 



II. THEORY 



A. Density in magnetic traps with finite depth 



The trapping potential due to a magnetic field B(r) with (local) minimum Bo is given by U(r) = mjgjfiB [B{r)—Bo] , 
where gj is the Lande g-factor and mj the magnetic quantum number of the state of the trapped atoms. We use 
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metastable helium in the 2 3 Si, mj = 1 state with gj w 2. The trapping potential of an IQ trap can be approximated 
by 



U[x, y, z) = y/a*(x* + V 2 ) + (U + /3z 2 ) 2 - U Q , (1) 

with minimum potential energy 

U Q = mjgj^ B Bo , (2) 

where the effect of gravity is neglected. 

According to Luiten et al. Q the atoms in a magnetic trap with a finite depth are described well by a Boltzmann 
energy distribution that is truncated at an energy equal to the trap depth. The phase-space distribution is assumed to 
be only a function of the energy of the atoms. At low temperatures this is the case if there are many elastic collisions 
between the atoms. At higher temperatures in the IQ trap, the atoms occupy a part of space large enough that 
higher-order terms of the potential break the axial symmetry of Eq. Q and the motion of the atoms becomes ergodic 
even without collisions. Below a summary is given of the equations needed to calculate the density distribution in a 
magnetic trap 0]. 

The (thermal) density distribution in an infinitely deep trap is 



noo(r) = n exp 



U(r) 



(3) 



k B T 

When only atoms with an energy smaller than e t are trapped, the density distribution becomes 

n(r)=P[§,K(r)]noo(r), (4) 

with incomplete gamma function P and 

^ (5, 
, U(r) > e t . 

In this truncated Boltzmann distribution T is still called the temperature, although strictly speaking a thermodynamic 
temperature is not defined for a nonequilibrium distribution. The central density is given by n(0) = P [|, fc(r)l no; 
parameter no is equal to the central density only in the limit of infinite trap depth. It is defined by 

n = N/V e , (6) 

with number of trapped atoms N and reference volume V e , which is equal to the effective volume N/n(0) in the limit 
of a deep trap. In general the reference volume is given by 

V e = A 3 C , (7) 

with thermal de Broglie wavelength 



/ 2tt^ 2 
y mkB-L 

and trap-dependent single-atom partition function £. For an IQ trap 

C = C^[P(4,r?) + |^P(3,n)], 

C = 6A 1Q (k B T)\ (9) 
A 1Q ^m 3 / 2 /(4VWa 2 h 3 ), 

with truncation parameter 



(10) 



The total internal energy of the trapped atoms is 



12P(5,n) + 6^P(4,n) 

E = -— k r f T Nk B T . (11) 

3P(4,n) + 2^P(3,n) 

Although no is not the central density and a thermodynamic temperature cannot be given for this nonequilibrium 
distribution, the (central) phase-space density (or degeneracy parameter) is still no^l 3 . 
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B. Evaporative cooling with constant cross section 



The behavior of a cloud of atoms during an rf sweep is simulated with the model of Luiten et al. |2| and Pinkse et 
al. 3] . The version of the model described here assumes that every atom with an energy greater than the trap depth 
e t leaves the trap (three-dimensional evaporation). Their trajectories should bring these atoms sufficiently fast to the 
exit area of the trap (one of the saddle points or, when an rf field is applied, positions where the rf field is resonant 
with the magnetic field), so they are removed without first colliding with another atom. Another assumption is that 
collisions occur in the s-wave regime. 

At the start, N atoms at temperature T are contained in a trap with depth e t determined by the magnetic field 
configuration. After a time step of negligible size the rf power is turned on: the rf frequency w r f, the time dependence 
of which is known beforehand, determines the truncation energy et = hiu r {— Uq. Then in small steps (for our situation 
10 ms) the loss of atoms and energy due to inelastic collisions, trap changes (spilling) and evaporation is determined, 
as discussed below. The temperature is found by solving numerically the equation for the total energy of the trapped 
atoms, Eq. 1|11[1 . If experimentally an extra temperature increase is found, for instance due to instability of the 
power supplies, this can be added separately. The calculation is stopped when either no atoms are left or when 
BEC is reached. In the last case the phase-space density is n A 3 > 53/2(1) 2.6, with polylogarithm function 

References y> y| do not give all loss rates in a suitable form, so first some intermediate equations are given and the 
energy density is derived. Trapped atoms have energy 

e (r, p) = U(r) + p 2 /2m . (12) 

The number of atoms with energy between e and e + de is p(e)f{e) de, with phase-space distribution 

«3. I £ 



/( e )=n yl d exp(-— J (13) 

and (in an IQ trap) an energy density of states 

P (e) = A 1Q (e 3 + 2Uoe 2 ). (14) 

The phase-space distribution can be given as a function of r and p via Eq. Ill2[l and is normalized so that the total 
number of trapped atoms is 

N '^II^ p)d ' ri ' p - (15) 

where the integration is done over the volume in phase space where e < £t- The density is given by 

n(r) = — - — ? / f(r,p)d 3 p 



^3 ^2m[e t -U(v)] , ^( r ) +p 2/ 2m 



exp 



(2-nhy Jo \ k B T 

= P[|,«(r)]noexp(-^). 



Anp 2 dp (16) 



Similarly, the energy density is 



-73 / e(r,p)f(r,p)d 3 p 
nn) J 



{p B TP[^K(r)]+U(r)P[l,K(r)]}n exp[- 



(17) 
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1. Inelastic collisions 

Inelastic collisions occur with background gas atoms, between pairs of trapped atoms and between three trapped 
atoms. Corresponding loss rates are 

N hgr = ~ J n(r) d 3 r = -N/t , (18) 
N 2h = -G Jn 2 (r)d 3 r, (19) 

A 3b = —L J n 3 (r)d 3 r, (20) 
respectively; corresponding energy loss rates are 

E hgI = - \ J e(r) d 3 r = -E/t , (21) 

E 2b = -G J e(r)n(r)d 3 r, (22) 

E 3h = -L J e(r)n 2 (r)d 3 r, (23) 

respectively. A dot denotes a derivative with respect to time. The integrals have to be calculated numerically for 
every time step. The constant r is the lifetime of the trap, determined experimentally. Fedichev et al. [T3. [T^ have 
calculated the other two constants for spin-polarized metastable helium. For T < 0.1 mK, G is only dependent on 
magnetic field. The maximum value is 6 x 10 -13 cm 3 /s at 750 G, but as the contribution of two-body collisions is 
more significant towards the end of the sweep when the cloud is concentrated at the center, the value for B < 10 G, 
which is G = 3 x 10 -14 cm 3 /s, can be used [l^. The three-body loss-rate constant is given by L = ll.6ha 4 /m, with 
scattering length a 13]. 



2. Spilling 



When the trap shape or depth is changed, the eigenstates of the trapping potential with highest energies can 
become unbound. Atoms in these states are spilled from the trap. This process does not depend on collisions; when 
the potential is changed only by lowering the truncation energy e t , spilling does not alter parameters T and Uq. 
After instantaneous lowering of the trap depth from e t to s' t , the change in the number of atoms due to spilling is 
AN S = — JJ p(e)f(s) de and the corresponding change in energy is AE e — — ep(e)f(e) de. With n A 3 — N/( 
and Eq. ©, integration yields 



ANf> = -N 



1 - 



3P(4,r/) + 2 <^P(3,r/) 



3P(4,r;) + 2^P(3,77) 



= — iV 



1 - 



AEff = -Nk B T- 



12 P(5, rj) + p( 4 , v ) - 12 P(5, V ') - 6^ P(4, 77') 



3P(4, ?7 ) + 2^P(3,7,) 



(24) 
(25) 



where rj' — e[/k^,T < 77; the values of N and T are those before the step takes place. 



3. Evaporation 

After an elastic collision of two trapped atoms with energy e < et, one atom may have an energy e > et and leave 
the trap. After this thermal escape, or evaporation, the average energy per atom has become smaller. The number 
and energy loss rates due to evaporation are 



^v = -l/ — n§<re-'7 CT , (26) 

V 7TTO 

E cv = N cv L + J k B T , (27) 
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with effective volumes for evaporation 

V cv = A*&, + Eli p (h v)] P(5, r?)} , (28) 

W„=A 3 £ 3 {(l + [v-EliHhV)] -P(M)} , (29) 

and elastic scattering cross section c The model has been developed for an energy-independent cross section a — 87ra 2 . 



C. Incorporation of an energy-dependent cross section 



The model as described in the previous section is restricted to the low-temperature limit, where a is energy- 
independent. In this section the theory will be adapted to include energy-dependence, although remaining in the 
s-wave regime. In our experiment the temperature at the start of the rf sweep is about 1 mK, where the dependence 
of a on the relative velocity of the colliding atoms cannot be neglected. For the derivation of atom loss rate iV cv 
and energy loss rate E cv with an energy-dependent cross section a, Luiten's derivation [14j is adapted. The elastic 
collision event rate r c , which is the number of collisions per second occurring in the whole cloud, is also derived. The 
average collision rate per atom is 2T C /N. This rate is not needed in the sweep simulation, but it is an interesting 
quantity, related to the atom loss rate N cv , which is the collision event rate under the condition that afterwards one 
of the atoms has enough energy to leave the trap. 



1. Collision event rate 



First the collision rate in the case of an (untruncated) Boltzmann velocity distribution is examined. The two-body 
elastic collision event rate is 

r c = [ |n(r) 2 (au r )d 3 r, (30) 



where rate coefficient (<jv t ) is the product of cross section a and relative velocity v T of the colliding pair of atoms, 
averaged over the thermal velocity distribution 

\ 3/2 / 2 

m \ 2 / UlVy 



^exp -77-^ )dvr. 



Assuming the rate coefficient is not dependent on position due to the magnetic field, the collision event rate can be 
written as 

r c = —^rn {<TV r ) (31) 



and the average collision rate per atom as (V2 e /V e )no(av T ), with 

n ( r ) j3„ 



V e = d'V = N/n , (32) 

J n 

V 2c =j(^±)\\. (33) 

For a Gaussian shaped density distribution Vi a jV c = 1/ v8, for the density distribution in an IQ trap with [To = (and 
infinite 77) V2 e /V e — 1/V32. The collision rate at position r is n(r)(av r ). Using the approximation a = 8ira 2 / \k 2 a 2 + 1) 
with thermal wave vector k = mv r /(2K), we get 

(av T ) = 8Tra 2 (v r ) [& - £^(0, &)] , (34) 



with shorthand notation £ c = ft 2 /(a^mk-oT) and average relative velocity (w r ) = Ay/k-gT/ (-Km). Figure^shows the 
quantity in square brackets and the temperature dependence of the rate coefficient in the case of metastable helium. 
In the low-temperature limit lim^oo^ — £ 2 e£r(0, £)] = 1. 
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FIG. 1: (a) Factor (av T ) / (8ira 2 (v r )) as a function of f, and (b) the rate coefficient as a function of temperature for metastable 
helium assuming four different values for the scattering length a. 



For a truncated Boltzmann distribution the collision event rate can be written as 

1 



2m(2irhy 



<r(\P2 -Pi\) \P2 ~Pi\f(r,pi)f(r,p 2 )d pi d p 2 d 



(35) 



The integration is over all (r,pi) with U(r) +pf/2m < et for both i = 1 and i = 2. A more convenient coordinate 
system uses average momentum P = (pi + P2)/2 and relative momentum q = P2 Pi, with angle between P and 



q and azimuthal angle With Jacobian d 3 pi d 3 p2 = d 3 Pd 3 g we get 



r, = 



2m(2nh) 

where the integration range is now given by 



J e -2U/k B T f {q)qe - q */4 m knT J ^PV^^T^^^^S 



(36) 



P 2 + -j + Pq\cos8\ < 2m[e t - U(rj\ 



(37) 



To eliminate the angle dependence, the integrand is multiplied by the fraction F of colliding pairs at position r for 
fixed P and q that are part of the truncated Boltzmann distribution: 



F 



1 

-in 



sin Odd d(p . 



where the integration over 9 is restricted by condition l|37l) . The result is 

Q 2 -P 2 -q 2 /A 



F = 



Pq 



(38) 



(39) 
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with Q — y/2m[et — U (r)] the maximum momentum an atom can have at position r without leaving the trap. This 
fraction has to be restricted to the physical range between and 1. For P > (Q 2 — q 2 /4) 1 / 2 at least one of the colliding 
atoms would have more momentum than is possible at position r and the fraction is zero. For < P < Q — q/2 all 
angles 6 are possible and the fraction is one. The integration over P is therefore divided into two parts: 



n 2 A 6 
2m{2irhf 



-2U/k B T 



2Q 



<T(q)qe~1 2 / imk * T 
rQ-9/2 



c -P 2 /mfc B r d 3p 



VQ 2 -? 2 /4 



0-8/2 



Fe -P 2 /mk B T d 3p d 3 qd 3 r ( 40 ) 



For an IQ trap [Eq. Q)] , the integration over position can be made one-dimensional Q using 

T[U{r)}d 3 r-. 4 " 



a 2 VPJo 



U(U + U Q )T{U)AU. 



(41) 



Introducing scaled variables y = q/ \frnk^T and «; = (e t — U)/k&T, and scaled constant 77 = et/k^T the collision 
event rate after integration over P becomes 



_ _ n 2 e- 2 *> 

J- r 



a 2 y/f3m 



(k B T) 2 



JO 



2 - 2cxp (y2ny - \y 2 ^j + y^yexp (2k- \y 2 ^j erf U/2k- \ 



dydn. (42) 



For an energy-dependent cross section a with hk = q/2 this equation has to be integrated numerically. For a 
constant cross section, after integration over y 7 the expression for r c contains terms 



f (77 - nfe lK erkVl^dn = T{x + l)c ln [P(as + 1 ,lrj) - P(ar + %,lrj)]/l x+1 



(43) 



and 



/ (77 - n) x ne K erfc V^d K = r(x + l)e r '{(r] + 1) [P(x + 2, 77) - P(x + |, 77)] 
./o 

- § P(x + §,77) - (as + 2) P(x + 3, 77) + (x + |) P(x + §,77)} , (44) 

with Eul er g amma function T(z) = r(z,0) and complementary error function crfc(z) = 1 — erf(z), which has the 
property 



e K erfc^=^ / 
t Jo 



Vt{t + K) 



dt. 



(45) 



The collision event rate for constant a becomes 

P c = &r$*y/*£A*C 00 {l + 8e-(r7 - 1 + §^77) 

7/ 4 + 477 3 -|-6?7 2 -6?7-27+|^(2?7 3 + 677 2 + 67;-3) j. (46) 



3 e 



When tr is energy- dependent, r c needs to be calculated numerically. However, in the limit of large 77 the velocity 
distribution becomes thermal and Eq. (f^ is expected to produce the same result as Eq. l|3T|) . For common values of 
77 and a cross section given by a = &na 2 / (k 2 a 2 + 1), Eqs. I|31|) and (|34|l can be used as an approximation: comparing 
exact numerical calculations of ^/(^"y^ 710 ) with the rate coefficient as given by Eq. I|34(l . at 77 = 10 deviations are 
smaller than 0.1% and at 77 = 5 smaller than 6% (for m < 100 u, \a\ < 100 nm, C/o/mb < 100 G, T < 5 mK). 



2. Evaporation rates 

The collision event rate r c is given by Eq. H36|) . with constraints to the integration range incorporated by a factor 
F in the integrand. The atom loss rate N ev is given by the same equation, but with an extra factor for the fraction 
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of collisions for which afterwards one of the atoms has enough energy to leave the trap. Since P and q are the same 
before and after collisions and all scattering angles 9' have equal probability in the s-wave regime, this factor is (1 — F). 
In the case of the collision event rate, the integration over P had to be divided into two parts [see Eq. l|40|l ] ■ Here 
only one part is left, as the integrand including F(l — F) with F — 1 is zero. The evaporation rate becomes 



a 2 y/P2m(2Trh) e J 



U(U + Uo)e~ 2U ^ T a(q)W^ q,i 



2Q 



y/Q 2 -q 2 /i 



-9/2 



F(l - F)4irP 2 e- p2 / mkBT dPdqdU , (47) 



and after integration over P: 



8a 2 ^/]3rn 



(k B T) s 



Jo 



u 

k B T 



<r(y) Vv-k[v - K + 

|-2[y(2 + y 2 - 8k) + 2\/2^(6 + y 2 - 8k)] exp (y/2Ky- ±y 2> 



+ 2 [8y + y/8K - y 2 (6 + y 2 - 8k)} - y/r [12 + Ay 2 + y 4 - 16(2 + jj 2 )k + 64k 2 



exp (2k - \y 2 ) [erf (J 2k - \y 2 ) - erf (7^- \y)\ \dydK. (48) 



For the energy loss rate E cv there is an extra factor in the integrand for the energy of the lost atom, 



1 



s 4 = U+-^[P 2 + + Pq\co$6' 
2m \ 4 

This means E ev is given by Eq. I|47|l with the factor (1 — F) replaced by 



G = T 



271 



e 4 sin 9' d8' dip' 



where the integration range is given by £4 > e t . After integration over 9' and if' we get 

(2P + q) 2 -4Q 2 ' 



G= (1-F) e, 
The energy loss rate after integration over P becomes 



16m 



(49) 



(50) 



(51) 



128a 2 



(k B T) 4 



77 pvSk 



Jo 



Up 



|256y + 2y / 8K~y 2 [60 + y 4 - 16y 2 (K - 2) + 64k(k - 1)] 

- 2{y 5 + 2y 4 V2k~ - 4?/ 3 (4k - 3) - 16y 2 \/2^ (2k - 1) 
+ y[68 + 32 K (2K-3)] +8V2^[15+ 16k(k- 1)]} exp (V2K~y - ±y 2 ) 
/tt (y 6 - 2y 4 (l2K -7)+ Ay 2 [21 + 8k(6k - 5)] + 8{15 - 4k[9 + Ak(Ak - 3)]}) 



exp(2K-i 2 / 2 )[erf[ v /2K-iy 2 J - erf [^-ij/JjjdzjdK. (52) 

When a is a constant, Eqs. I|48|) and H52fl reduce after integration to Eqs. (|26|l and \21\ . respectively. When a is 
energy-dependent, the expressions for 7V CV and E cv have to be integrated numerically. An effective cross section a e g 
can be introduced, which will give the correct value of N ev if a = a c g in Eq. 126|) . In the case that the cross section 
is given by a — 8na 2 / (k 2 a 2 + 1), the effective cross section can be approximated by a variation of Eq. 



a cS « 8™ 2 [Ccv - &e«-T(0, £ ev )] , 



(53) 
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FIG. 2: Effective cross section as a function of temperature for m = 4 u, a — 10 nm, Bo = 0.46 G and four values of r) (solid 
curves), together with the approximation (dashed curves). The dotted line is the low-temperature limit. 



where 



Scv 



3h 2 



r] mk^Ta 2 Etma? 



(54) 



In this approximation, found by trial and error, the effective cross section is only dependent on temperature via the 
trap depth. Figure [3 shows a c g and its approximation as a function of temperature in our situation for a — 10 nm 
and several values of 77. Apparently a c g(T) does not approach an asymptote at high 77. In general (m < 100 u, 
\a\ < 100 nm, C/q/mb < 100 G, T < 5 mK, 2 < 77 < 25) the absolute difference between a e g and its approximation is 



less than 3% of the low-temperature limit 



shown as a dotted line in Fig. [3 The energy loss rate can still be 



calculated with Eq. 
2%. 



I|27|l . even though the effective volumes are not correct: the relative error in E cv /N cv is less than 



III. EXPERIMENT 

We have performed evaporative cooling experiments in a setup which is described in detail by Herschbach et al. |(| . 
In short, an atomic beam from a DC discharge source of metastable helium is collimated, deflected and slowed in 
a traditional Zeeman slower. The atoms are trapped in a magneto-optical trap (MOT), cooled in optical molasses, 
spin-polarized by optical pumping, again trapped and finally compressed in a cloverleaf magnetic trap, an example 
of an IQ trap. The magnetic field configuration is determined by a/(2/ie) = 69 G/cm, /3/(2/jb) = 13.9 G/cm 2 , 
Bq = 0.5 G, and the trap depth is 47 G (corresponding to 6.4 mK). The harmonic part is characterized by trap 
frequencies lo p /(2tt) — 853 Hz and ui z /{2n) — 44 Hz. After compression a cloud of ~10 9 atoms at a temperature of 
~1 mK and a central density of ~10 10 cm~ 3 is obtained. Temperature T and a relative measure of the number of 
atoms are determined by time-of-flight (TOF) measurements with a microchannel plate (MCP) detector positioned 
18 cm from the trap center, whereas calibration of the number of atoms N is performed by absorption imaging in the 
MOT. As the radial confinement is much stronger than the axial confinement, the cloud is elongated: the theoretical 
density while the atoms are still in the trap has a FWHM of 0.2 cm horizontally and 1.4 cm vertically. A contour 
plot of the density, integrated in one horizontal dimension, is shown in Fig. |3Ja). Absorption images are taken 1 ms 
after the trap has been switched off. To simulate the ballistic expansion during this time, the theoretical density 
distribution is convolved numerically with a Maxwell-Boltzmann velocity distribution before the integration in one 
dimension, using a temperature of 1.2 mK determined from a TOF measurement [Fig.|3{b)]. No fitting is required and 
the shape agrees well with the experimental column density [Fig. Etc)]. In contrast, a Gaussian fit to the experimental 
column density deviates significantly [Fig. Old)]. At lower temperatures the cloud becomes smaller and a larger part 
is confined to the harmonic region of the trap. Therefore the density distribution approaches a Gaussian shape during 
evaporative cooling. However, only below 17 /iK will the r.m.s. radii calculated with a harmonic approximation of 
the trap deviate less than 10% from the exact solution. 

For evaporative cooling, the rf frequency starts at 120 MHz and decreases exponentially with a time constant of 5 s, 
close to the optimum found by maximizing the elastic collision rate per atom after a sweep to 5 MHz. The progress 
of the temperature and number of atoms is followed both in absorption imaging and with the MCP detector. A series 
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FIG. 3: The density distribution integrated in the x-direction. From left to right: the theoretical distribution (a) in the 
cloverleaf trap and (b) after 1 ms expansion, (c) the density from the experimental absorption image (the negative logarithm 
of the transmittance), and (d) a Gaussian fit to the experimental density distribution. Contours at 10%, 30%, 50%, 70% and 
90% of the maximum integrated density. 
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FIG. 4: Collection of absorption images taken 1 ms after the cloverleaf trap is switched off, at various stages during an rf sweep. 
From left to right: just before rf is turned on, at 100 MHz, 65 MHz, 35 MHz, 15 MHz, 10 MHz, and 5 MHz. 



of images taken after 1 ms expansion is shown in Fig.0] To emphasize the change in cloud size, maximum absorption 
within each image is rendered black. The central absorption varies from 62% in the leftmost image to 14% in the 
rightmost image. The smallest cloud contains 9(4) x 10 6 atoms at a temperature of 23(3) /iK. Above 0.5 mK the radii 
are larger than expected, especially in the horizontal plane, with a density that is half of the predicted value. This 
may be explained by a time lag between a change of trap depth and getting closer to steady state at a smaller cloud 
size: the sweep has to be executed rather fast, due to a limited trap lifetime of 12.5 s. For theoretical calculations the 
actual size during the sweep is required. A first order correction of the modeled density distribution is to multiply 
trap parameters a and j3 in Eq. Q with factors that are kept constant during the sweep. These are f a = 0.74(9) and 
fp = 0.89(5), respectively, where the uncertainty includes the variation in the measured radii and the temperature 
uncertainty. 

Figure shows the temperature and number of atoms as a function of time since the start of the rf sweep. Each 
point is the result of a TOF measurement after the sweep is interrupted by turning off the rf power and the trap at the 
same time. Error bars are given as a gray band, including a 20% standard deviation in the calibration of the number 
of atoms according to the MCP and including the uncertainty in the temperature due to magnetic field gradients 
between the trap and the MCP ^(|. The indirect effect of the temperature error on the number of atoms is also 
taken into account. The line is a theoretical calculation, assuming the scattering length is 10 nm. The trap lifetime 
is measured to be 12.5 s, the measured heating of 0.7 /iK/s is neglected, and the trap minimum is 0.5 G. Correction 
factors f a and fp are included, so a/(2t*B) = 51 G/cm and /3/(2/xb) = 12.4 G/cm 2 . Due to an unexpected feature of 
the rf generator, below 7 MHz the rf power drops, dwindling to zero below 4 MHz. This means the rf is effectively 
turned off after 17 s, the temperature stays constant and the number of atoms goes down only due to background 
collisions. 

The five loss rates used in the model are plotted as a function of time in Fig. H3 for the theoretical sweep with 
a = 10 nm. The contribution of two- and threediody collisions can be neglected up to the moment BEC is reached, 
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FIG. 5: (a) Temperature and (b) number of atoms as a function of rf time during a sweep. The gray bands show the possible 
values due to uncertainty in the detector calibration and in the temperature. The curve is a simulation with a — 10 nm. 

theoretically after 22 s. At the start of the sweep the dominant loss process is spilling, because the frequency is 
decreased relatively fast. If the sweep would take longer, spilling losses would be reduced in favor of evaporation 
losses, but this would only be advantageous if the trap lifetime were larger. Figure EJb) shows that an exponentially 
decreasing rf frequency goes unnecessarily fast towards the end: the spilling contribution increases in the last seconds. 

To find the scattering length best describing the data, one can simulate sweeps for several values of a starting from 
the temperature and number of atoms in the beginning. However, deviations early on influence the behavior of the 
curve at later times. For instance, turning on the rf is not well modeled, as the experimental temperature after 1 s is 
lower than expected [see Fig. 1ST a)] . In addition, the temperature as a function of time exhibits a slope change after 
lis, which is not seen in simulations with any value of the scattering length. Therefore, small sweeps are simulated 
for each interval between the data points. Starting from the experimental situation at one point with a constant a 
in Eq. J2HJ), the temperature at the end of the sweep is calculated. The values of a for which the end temperature 
according to the model corresponds to the experimental temperature are given in Fig. EJa) . The curve is the effective 
cross section according to Eq. (|53|l for a = 10 nm. The gray band indicates the error due to the uncertainty in 
the MCP calibration and the temperature (as in Fig. |3J) , as well as the uncertainty in the correction factor for the 
theoretical cloud size. The scattering length determined with the experimental effective cross section and Eq. 15311 is 
given in Fig. Efb). The point in the first interval is dropped, because the effective cross section was unusually high. 
This in turn is caused by the unexpectedly large drop in temperature after turning on the rf power, which makes 
evaporation look very effective. In the last part of the sweep the uncertainty in the scattering length becomes large 
due to the increasing uncertainty in the temperature determination. However, the possible error in a in the first half 
of the sweep may be larger than given, because the evaporation process seems not to behave entirely according to the 
model. Also, at the start we have to rely more on the model anyway, as the scattering length is determined only from 
evaporation, which then constitutes just a quarter of the loss rate [Fig. Efb)]. Taking this into account a scattering 
length a = 10(5) nm can be considered as consistent with the experimental data. 
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FIG. 6: (a) Absolute and (b) relative loss rates due to background collisions M> gr (solid curve), evaporation N cv (long-dashed 
curve), spilling Ng (dotted curve), two-body collisions A^b (dash-dotted curve) and three-body collisions A^b (small-dashed 
curve), as a function of sweep time. In (b) the three loss rates due to bad inelastic collisions are summed. This is according to 
a simulation with a = 10 nm. 

IV. SUMMARY AND CONCLUSIONS 

To summarize, we have developed a model that allows simulation of evaporative cooling experiments in situations 
where the elastic scattering cross section depends on collision energy. The model is used to simulate the number of 
atoms and temperature in an rf sweep for metastable helium. From a comparison with the experiment an s-wave 
scattering length \a\ = 10(5) nm is extracted. This value is a factor of two smaller but within the experimental 
accuracy of previous experiments; the BEC experiment of the Orsay group yields a = +20(10) nm 0,^3' ^ ne result 
of the ENS group is a — +16(8) nm 0,13. Later experiments of the ENS group are more consistent with a lower 
rather than a higher value 0. The scattering length can also be determined theoretically, because the 5 E+ potential 
in which spin-polarized atoms collide has been calculated. From the Starck and Meyer potential ||| a scattering length 
a = +8.3 nm is deduced. More recently, Gadea, Leininger and Dickinson ^(J calculated the short-range part of this 
potential more accurately and combining their potential with the accurately determined long-range potential of Yan 
and Babb 19] they determine a = +15.4 nm. However, the value of the scattering length is very sensitive to the way 
the short-range and long-range potential are connected. An improved method of combining the two potentials results 
in a = +12.2^4 1 nm Our experimental result is in good agreement with this value. 
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FIG. 7: (a) Effective cross section and (b) scattering length as a function of sweep time. The curve is a simulation with 
a — 10 nm, the gray bands indicate the error due to uncertainty in the number of atoms, temperature and cloud size. 
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